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To study the properties of frustrated two-dimensional Josephson junction arrays near incommen- 
surability, we examine the current-voltage characteristics of a square proximity-coupled Josephson 
junction array at a sequence of frustrations / = 3/8, 8/21, 0.382 (« (3 — \/5)/2), 2/5, and 5/12. 
Detailed scaling analyses of the current-voltage characteristics reveal approximately universal scal- 
ing behaviors for / — 3/8, 8/21, 0.382, and 2/5. The approximately universal scaling behaviors and 
high superconducting transition temperatures indicate that both the nature of the superconducting 
transition and the vortex configuration near the transition at the high-order rational frustrations 
/ = 3/8, 8/21, and 0.382 are similar to those at the nearby simple frustration / = 2/5. This finding 
suggests that the behaviors of Josephson junction arrays in the wide range of frustrations might be 
understood from those of a few simple rational frustrations. 
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For Josephson junction arrays (JJA's), frustration can 
be introduced by applying an external magnetic field. In 
the presence of frustration, a finite density of vortices is 
induced in the array. For a two-dimensional (2D) JJA, 
the mean number of vortices per plaquette induced by 
a uniform external field is identical with the frustration 
index / defined as the fraction of a flux quantum per 
plaquette. The competition between two different length 
scales, the mean separation of vortices and the period of 
the underlying pinning potential of the array, leads to 
various structures of vortex lattices at low temperatures. 
At a simple rational /, the ground state is a commen- 
surate pinned vortex latticed For example, vortices in a 
square array at / = 1/2, 1/3, or 2/5 have a staircase 
state with a quasi-one-dimensional structure in a diago- 
nal direction.'^ For /= 1/3 or 2/5, the high-temperature 
isotropic vortex liquid freezes into the pinned vortex lat- 
tice with the staircase configuration through a finite- 
temperature superconducting transition. In the weak 
magnetic- field limit, the superconducting vortex lattice 
becomes triangular.'^ In a high-order rational or an ir- 
rational field, however, both the low-temperature vor- 
tex configuration and the nature of the superconducting 
transition remain controversial. Some years ago, Halsey^ 
suggested, based on Monte Carlo (MC) studies of the XY 
model, that in the limit of the irrational /* = (3 — \/5)/2, 
there may exist a glass transition to a superconducting 
disordered vortex state at some finite temperature Tg « 
0.25 J//cb for a square array, in which J is the junction 
coupling strength. Other argumentSf^*^ however, sug- 
gest that the glass transition appears at zero temperature 
and that the superconducting transitions at rational /'s 
near /* should occur at some finite temperatures which 
decrease to zero as / approaches the irrational value of 
/*. Some of the recent MC simulationSj^ on the other 
hand, showed that at most /'s in the range 1/3 </ < 
1/2 except a few simple rational /'s, vortices in a square 
JJA undergo two separate transitions at finite temper- 
atures: a sharp first-order phase transition at temper- 
ature To = (0.15— 0.2) J/fcs to a resistive quasiordered 



state and a pinning transition at a lower temperature to 
a superconducting ordered state with the vortex configu- 
ration consisting of periodic arrangements of either stair- 
case patterns of the simple rational /'s or hole defects on 
the f — 1/2 configuration. A single finite-temperature 
phase transition from a vortex-liquid phase to a supercon- 
ducting ordered phase consisting of domains separated 
by parallel domain walls has also been suggested by a 
different numerical work for a sequence of high-order ra- 
tional /'s approaching the irrational f *.^ The nature of 
the superconducting transition and the vortex structure 
of the superconducting ground state near incommensu- 
rability obtained by numerical efforts have been found^ 
to strongly depend upon the imposed boundary condition 
and choice of dynamics. 

There have been very few experimental studies about 
phase transitions near incommensurability. They agree 
on a single finite-temperature phase transition to a su- 
perconducting state which can be understood as a pinned 
vortex phase with long-range phase coherence jiSiii How- 
ever, since the experiments were performed on super- 
conducting wire networks where phase fluctuations were 
much weaker than in JJA's, the superconducting behav- 
ior is suspected by some to be dominated by a mean- 
field transition!^ Besides, the experiments provided only 
limited information about the vortex structure of the su- 
perconducting state. In addition to stronger phase fluc- 
tuations, a proximity-coupled JJA has a broader critical 
region than a superconducting wire network, which en- 
ables one to investigate the superconducting transition 
with higher precision. For a JJA, it is also possible to 
measure the single-junction critical current ic and deter- 
mine the phase-transition temperature in units of J/ks 
(= hic/2ekB)- Knowing the phase-transition tempera- 
ture in units of J/ks, one can make a direct compari- 
son between experimental results and the numerical ex- 
pectations. In this paper, we present an experimental 
investigation of the superconducting phase transition of 
a proximity-coupled JJA near incommensurability, that 
is, at high-order rational frustrations near the irrational 
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FIG. 1: Some of IV curves for four different frustrations: (a) 
/ = 3/8 at r = 5.150 to 5.900 K, (b) / = 0.382 at T = 5.200 
to 5.900 K, (c) / = 2/5 at T = 5.150 to 5.876 K, and (d) / = 
5/12 at T = 5.170 to 5.878 K. The dashed lines are drawn 
to show the power law {V ~ 7^"'"^) behavior at the critical 
temperature. 
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FIG. 2: Current dependences of the slope of the IV curves in 
Fig.0 



/* = (3 — •\/5)/2. We examine the superconducting scal- 
ing behavior of the current-voltage {IV) characteristics 
of a square array at / 3/8, 8/21, 0.382 (w (3- V5)/2), 
2/5, and 5/12. Detailed scaling analyses reveal approx- 
imately universal scaling behaviors for / = 3/8, 8/21, 
0.382, and 2/5. The superconducting transition temper- 
atures are found as high as (0.19 — Q.22)J /ks- We dis- 
cuss implications of the universal scaling behavior and 
the high superconducting transition temperatures in con- 
junction with the nature of the superconducting transi- 
tion and the vortex configuration of the superconducting 
state near incommensurability. 

The experiments were performed on a square array 
of 200x1000 Nb/Cu/Nb Josephson junctions described 
in Ref. 12. The array completed the zero- field super- 
conducting transition at T 6.24 K with a transition 
width « 0.4 K. The appearance of many resistance min- 
ima in the magnetoresistance curve confirms the good 
uniformity of the sample. The frustration could be pre- 
cisely adjusted by the use of the magnetoresistance curve 
of the sample showing sharp resistance minima at frac- 
tional f'sm^ The standard four-probe method utilizing a 
transformer-coupled lock-in voltmeter and a square- wave 
current at 23 Hz was used for the IV characteristics mea- 
surements. The single-junction critical current ic and the 
junction coupling strength J{= Uic/^e) at high tempera- 
tures were determined by extrapolating the ic vs T data 
at low temperatures by the use of de Gennes expressions^ 
in the dirty limit. The ic at low temperatures was ob- 
tained from the / vs dV/dl curves. 

The IV characteristics of the sample for / = 3/8, 
0.382, 2/5, and 5/12 are shown in Fig. ^ Relatively 
narrow voltage ranges of the IV data are related to the 



small dynamical critical exponent for frustrated JJA's. 
As can be seen from the current dependences of the slope 
of the IV curves, shown in Fig. |21 the IV curves at high 
temperatures for all the /'s are Ohmic at low currents 
and bend upward at high currents. At low temperatures, 
the IV curves arc bent downward. The low-temperature 
curves can be fitted into a form V ~ /exp[— (/t//)'^] 
with fi = 0.6—1. The activated character of the low- 
temperature curves proposes a superconducting state as 
the low-temperature state. The variation of qualita- 
tive nature of the IV curves suggests that for all the 
/'s, a phase transition from a high-temperature resis- 
tive state to a low-temperature superconducting state 
occurs at the temperature where a straight IV curve 
appears. The gradual drop of the resistance upon ap- 
proaching the transition is indicative of a continuous 
phase transition. For a continuous superconducting tran- 
sition in 2D, IV curves are expected to collapse onto two 
curves with / and V scaled by the general scaling form 
V/I\T - Tcl^" = £±{I/T\T - Tc]"), where z is the dy- 
namical critical exponent, the correlation-length crit- 
ical exponent, and £± the scaling functions above and 
below the transition temperature Tci^ This scaling form 
becomes a simple power-law IV relation V ^ j^+i at 
T ^ Tc and V/I - (T - TcY" in the low / hmit at 
T > Tc- Figure 121 shows the IV curves scaled with the 
scaling form. For all four /'s, we find the IV data ex- 
hibit good scaling behaviors. A good scaling behavior is 
found for / = 8/21 as well. The scaling parameters 
V, and z for / = 8/21 are practically identical with those 
for / — 0.382. In the scaling analyses, we made use of 
the approximate values of Tc, z, and v derived from the 
straight IV curves in Fig. ^ f^nd the temperature vs re- 
sistance curves fitted into the relation R ^ {T — Tc)^'^. 
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FIG. 3; Scaling plots of the IV curves. Each plot contains 
IV curves at 18—21 different temperatures. The values of Tc, 
z, and v used to scale the data are shown in the insets. 



The values of Tc, v, and z in the insets of Fig. O are the 
refined ones in the scaling processes. The error bars, that 
is, the arbitrarinesses in determining the scahng param- 
eters are 0.03 K for Tc, 0.05 for z, and 0.1 for v. Distinct 
concavities of the IV curves significantly limit the arbi- 
trariness in determining the parameters. The IV data 
satisfy the recently proposed criterioniS. to determine if 
the data collapse is valid, as appeared in Fig. |21 The 
effect of finite-size-induced free vortices, which may al- 
ter the scaling behavior considerably in the absence of 
a magnetic field'''^, is not significant for our sample ex- 
posed to a magnetic field. We thus conclude that the 
scaling plots confirm for all the frustrations investigated 
a superconducting transition at the temperature where a 
straight IV curve appears. 

The Tc's for the /'s investigated correspond to (0.19- 
0.22) J//cb, which is somewhat too high to be related with 
the vortex-pinning transition, the lower-temperature su- 
perconducting transition of two separate phase transi- 
tions observed in some MC simulationsi^ The vortex- 
pinning transition was expected to appear at Tc < 
QAJ/kB for / = 3/8 and 0.382. The weak / depen- 
dence of Tc does not seem to be compatible either with 
the argumen1i& that Tc should decrease to zero as / ap- 
proaches the irrational /*. The critical exponents v and 
z for / = 3/8, 8/21, and 0.382 are identical within exper- 
imental errors with those of / = 2/5. In Fig.01 we find 
that not only are the critical exponents similar for four 
different frustrations but the scaling functions are also. 
However, the IV data for / = 5/12 which is only ~ 0.017 
distant from 2/5 fail to collapse onto the same curves in 
Fig.H For / = 1/3, it was foundi^ that v = 1.6 and z = 
0.60, significantly different from those of / = 2/5. Such 
approximately universal scaling behaviors for / = 3/8, 
8/21, 0.382, and 2/5 indicate that the superconducting 
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FIG. 4: Scaling plot of the IV curves for / = 3/8, 8/21, 0.382, 
and 2/5 with the same z (= 0.76) and v (= 1.84). The IV 
data for four different frustrations collapse onto the nearly 
same curves. 



transitions for / ~ 3/8, 8/21, and 0.382 are similar in 
nature to that for the nearby simple frustration / — 2/5. 
For / = 2/5, the superconducting transition has been un- 
derstood in terms of the freezing of a vortex liquid into a 
pinned ordered vortex solid. ^^'■^"'^^ The approximate uni- 
versality also implies that the vortex configurations near 
Tc for / = 3/8, 8/21, 0.382, and 2/5 are quite similar 
to each other. Adsorbed atoms on a periodic substrate 
which bear a strong resemblance to vortices in a periodic 
array of Josephson junctions may have a ground state 
near incommensuracy consisting of domains of nearby 
commensurate phase (s) separated by domain walls.- ^ If 
the JJA has domains of the / = 2/5 phase separated 
by quasiperiodic or commensurate domain walls as the 
ground states for / — 3/8, 8/21, and 0.382 analogous 
to the adsorbate systems near incommensurability, the 
vortex configurations for / = 3/8, 8/21, 0.382, and 2/5 
may look similar in the vicinity of transition where many 
extra domain walls get excited. Extending our finding 
at the /'s around 2/5, one might expect that the behav- 
iors of vortices in the wide range of frustrations could be 
understood from those of a few simple rational frustra- 
tions at which the system has distinct local minima in 
energy. The IV characteristics for / — 5/12 might then 
be understood from those of possibly f — 3/7. 

As discussed above, our experimental finding is in- 
consistent with the numerical results indicating either 
a glasslike transition or two separate transitions at the 
irrational /* or its nearby high-order rational /'s. In- 
stead, it seems to be compatible to some extent with 
the simulation results^ indicating for a sequence of /'s 
approaching the irrational /* a single finite-temperature 
ordering transition of a vortex liquid into a striped lat- 
tice phase with domains of the / — 8/21 staircase state 
separated by parallel walls. Denniston and Tang^ argued 
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that the conflicting simulation results near the irrational 
/* are due to the impositions of different boundary con- 
ditions. An imposition of an improper boundary condi- 
tion incompatible with incommensurate or long-period 
commensurate phases may produce unrealistic simula- 
tion results near incommensurability. Another possible 
cause for the inconsistencies between the experiments 
and the simulations is the disorder effects. A real ar- 
ray always contains some amount of quenched disorder, 
for instance, the inevitable variation of the junction cou- 
pling strength. Even a small amount of quenched disor- 
der has been foundiiiiSiSiiSiSi to have a significant effect 
on phase transitions of JJA's at simple rational frustra- 
tions. Therefore, for the present, we may not totally 
exclude the disorder effects as a possible cause for the 



inconsistencies. 

In conclusion, approximately universal scaling behav- 
iors of the IV characteristics are found for a square JJA 
at / = 3/8, 8/21, 0.382 (3- V5)/2], and 2/5. The ap- 
proximately universal scaling behaviors and high super- 
conducting transition temperatures indicate that both 
the nature of the superconducting transition and the vor- 
tex configuration near the transition at the high-order 
rational frustrations / = 3/8, 8/21, and 0.382 are similar 
to those at the nearby simple frustration / = 2/5. This 
finding suggests that the behaviors of Josephson junction 
arrays in the wide range of frustrations might be under- 
stood from those of a few simple rational frustrations. 
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